Frequency Bin Entangled Photons 
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A monochromatic laser pumping a parametric down conversion crystal generates frequency entan- 
gled photon pairs. We study this experimentally by addressing such frequency entangled photons 
at telecommunication wavelengths (around 1550 nm) with fiber optics components such as electro- 
optic phase modulators and narrow band frequency filters. The theory underlying our approach is 
developed by introducing the notion oi frequency bin entanglement. Our results show that the phase 
modulators address coherently up to eleven frequency bins, leading to an interference pattern which 
can violate a Bell inequality adapted to our setup by more than five standard deviations. 

PACS numbers: 42.50. Dv, 03.67.Bg, 03.65. Ud 



I. INTRODUCTION 



Entanglement is one of the most fascinating aspects 
of quantum mechanics, used both for fundamental tests 
of physical principles and for applications such as Quan- 
tum Key Distribution (QKD). Many different kinds of 
photonic entanglement have been produced, including 
entanglement in polarization [l|, Q , momentum Q , an- 
gular momentum Q and time-energy. Investigation of 
the latter degree of freedom has been mainly inspired 
by two photon bunching experiments first carried out by 
Ou and Mandel [sj, see @, 0]; and by Franson's proposal 
Q for addressing the entanglement in the time domain, 
see [ol-fT^. Photons entangled simultaneously in both 
time-energy and other degrees of freedom have also been 
studied [3, [iBl ■ Time-energy entanglement can also be 
viewed as frequency entanglement, as demonstrated in 
recent works [ij, [l^. Here we show how to address 
time-energy entangled photons directly in the frequency 
domain. This is realized in optical fibers by using com- 
mercially available telecommunication components. 

Before presenting our approach, it may be useful to 
recall Franson's proposal Q which is based on three key 
ideas. First, a continuous pump laser produces time en- 
tangled photon pairs: the emission time of each photon is 
uncertain, but both photons are emitted simultaneously. 
Second, one uses measurements that resolve the arrival 
time of the photons. This leads to the concept of time 
bin: two photons whose arrival time cannot be distin- 
guished by the detectors belong to the same time bin. 
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Third, different time bins are made to interfere by using 
unbalanced Mach-Zehnder interferometers. These ideas 
provide a powerful platform to investigate quantum en- 
tanglement, yielding seminal works such as long distance 
violation of Bell inequalities [llj and entanglement based 
QKD [H. 

Here we also use the time-energy degree of freedom, 
but the way it is addressed is very different. However at 
the conceptual level there is an instructive parallel be- 
tween our approach and that of Franson. First, a narrow 
band pump laser produces frequency entangled photon 
pairs: the frequency of each photon is uncertain, but the 
sum of the frequencies is well defined. Second, our detec- 
tors are preceded by narrow band filters that resolve the 
frequency of the detected photons. This leads to the con- 
cept of frequency bin: two photons whose frequency is so 
close that they cannot be distinguished by the filters are 
said to lie in the same frequency bin. Third, different fre- 
quency bins are made to interfere by using electro-optic 
phase modulators. 

Our work is inspired by, or related to, earlier proposals 
for manipulating qubits in the frequency domain |19l425l | . 
Our experimental techniques follow closely those of QKD 
systems in which the quantum information is encoded in 
frequency sidebands of an attenuated coherent state [l9l - 
[2l| . Such systems allow efficient transmission of quan- 
tum information at telecommunication wavelengths. The 
main advantage of this method for encoding and carrying 
out transformations on optical qubits is that one does not 
need to stabilize paths in optical interferometers. Rather 
one must only lock the local Radio Frequency (RF) oscil- 
lators used by Alice and Bob, which is much easier. Fur- 
thermore information encoded in sidebands is unaffected 
by birefringence in the optical fiber used for transmission. 
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Recent improvements to these experiments have included 
dispersion compensation and l ong distance synchroniza- 
tion of the sender and receiver [221 ' so that this approach 
constitutes now the only commercial alternative to time 
bin based QKD. 

The architecture reported in [l^, [23| was dedicated to 
QKD using faint laser pulses, but it is inefficient when 
single photons are used because weak modulation am- 
plitudes are required. To overcome this limitation, an 
alternative method was proposed |21| in which informa- 



tion is encoded both in the amplitude and relative phase 
of three frequency bands generated by electro-optic phase 
modulators. This second approach is attractive because 
in principle the phase modulators need not attenuate 
the signal, because there is no need for a strong refer- 
ence pulse, and because the phase modulators can ad- 
dress many frequency sidebands simultaneously. Here 
we transpose, with appropriate modifications, the setup 
of [111 to the entangled photon case. 

In the following we first describe our experiment and 
give the principle of our method. Then we present our ex- 
perimental results, and demonstrate that the phase mod- 
ulators can coherently address up to eleven frequency 
bins. Note that high dimensional entanglement has 
been studied in a number of earlier experiments, see e.g. 
[H [H, [H-ll^. We finally show that the two photon 
interference pattern we obtain can in principle violate a 
Bell inequality adapted to our setup by more than five 
standard deviations. 



II. EXPERIMENTAL SETUP, THEORETICAL 
DESCRIPTION, AND RESULTS 

Our all fiber experiment is schematized in Fig. [T] 
Photon pairs are generated by parametric down con- 
version in a 3 cm long Periodically Poled Lithium Nio- 
bate (PPLN) waveguide (HC Photonics). PPLN waveg- 
uides have emerged as the preferred photon pair source 
at telecommunication wavelengths because of their ex- 
tremely high spectral brightness [3l|, H^. The narrow 
band pump laser (Sacher Lasertechnik, Ap = 773.865 nm, 
Pp « 6 mW) is removed with a drop filter F insuring 
125 dB isolation. The identically polarized photon pairs, 
distributed around Ao = 1547.73 nm = 2Trc /luq, are sepa- 
rated with a 3 dB coupler. Interesting cases occur when 
the photon pair is split: one photon is sent to Alice (A) 
and the other to Bob (B). At the output of the coupler, 
the photons pass through electro-optic Phase Modulators 
PMa,b (EOSPACE, 25 GHz bandwidth, 2.5 dB loss, half- 
wave voltage Vtt w 2.9 V), whose active axis are aligned 
with the linear polarization of the photons (preserved 
thanks to polarization maintaining fiber components), 
and to which are applied sinusoidally varying voltages 
at frequency ri/27r = 12.5 GHz, with amplitudes Va.b 
and phases a, /3 which can be controlled. The induced 
time dependent optical phases 4'Ait) = acos{Qt — a) 
and (pBit) = bcos{nt — /3), where a — ttVa/K- and 



b = ttVb/Vtt, lead to the unitary transformations 

\uj)a ^\uj +pn)AUp{a,a) , (1) 



Mb 



E 



LJ + qn)BUg{b,(3), 



(2) 



where subscripts A and B refer to Alice and Bob's 
photons, Up{a,a) = ,/j,(a)e*P("-'^/2), Ug{b,P) = 
Jg(6)e"^*^'^~'^/^\ and Jp^q is the p, gth-order Bessel func- 
tion of the first kind. The range of values of a, b which 
were experimentally accessible are limited to {0, 2.74} 
due to the finite power of the RF generator used. The 
photons are then sent through narrow band filters Fa,b 
which consist of a Fiber Bragg Grating (FBG) preceded 
by a circulator. Losses are 0.2 dB for the FBGs and 
0.8 dB for the circulators (round-trip). The spectral char- 
acteristics of the FBGs are: Full Width at Half Maximum 
(FWHM) sa 3 GHz, and more than 30 dB isolation at 
6.25 GHz. see Fig. [l]inset (a). Alice's filter is kept fixed 
on angular frequency uja = ^o- It is athermally pack- 
aged to reduce central wavelength deviation to 1 pm/K. 
The temperature of Bob's filter is controlled by a Peltier 
module, which allows continuous tuning of the reflected 
frequency ujb over a 1 nm range around ujq . The use 
of such narrow band filters together with a spectrally 
bright PPLN source of entangled photons has been re- 
ported previously in the context of four-photon exper- 
iments HI, 111]. Finally the photons are detected by 
two Avalanche Photo-Diodes APDa,b operated in gated 
mode (id Quantique, efficiency 15%, dark count rates 
3.5 • 10~^/ns and 8.0 • 10~^/ns) and a time to digital con- 
verter performs a coincidence measurement. The maxi- 
mum coincidence rate was approximately 10 Hz, which is 
consistent with the earlier work of [3^ . [33| . 

A parametric down conversion source pumped by a 
monochromatic beam produces an entangled state which 
we can idealize as 



(3) 



The total energy of the photon pair is well defined, but 
the energy of each photon is uncertain. For simplicity of 
notation we have not normalized Eq. ([3]). This does not 
affect our predictions as we are interested in the ratios 
of the probabilities of finding photon A at one frequency 
and photon B at another frequency for different settings 
a, 6, a, /3 of the phase modulators. For a discussion of how 
to normalize Eq. ^ so as to describe a rate of photon 
pair production, see the appendix. 

Note that taking the Fourier transform of Eq. ([3]) 
would yield a description of the state in terms of time 
entanglement: the arrival time of each photon is uncer- 
tain, but the difference between the arrival time of Alice 
and Bob's photon is well defined. The approximations 
leading to Eq. ^ consist in neglecting the finite pump 
bandwidth (which is approximately 2 MHz) and the finite 
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signal and idler photons bandwidths (which are approx- 
imately 5THz). This is legitimate as they are respec- 
tively much smaller than the bandwidths of the filters 
Fa,b (which are approximately 3 GHz, see Fig. [T] inset 
(a)) and much larger than the bandwidth sampled by the 
phase modulators (which is approximately f 25 GHz, see 
Fig. [U inset (b)). 

According to Eqs (|f I2p . the phase modulators realize 
interferences between photons whose frequencies are sep- 
arated by integer multiples of fi. They thus play the same 
conceptual role as the Mach-Zehndcr interferometers in 
Franson's scheme which realize interferences between dif- 
ferent time bins. Using Eqs (|1I2[) one can readily compute 
how the entangled state Eq. ^ is affected by the phase 
modulators: 

I*) ^ / duj'^\uJo + uj')A\i^o-uj' + dn)BCd{a,b,a,l3), 

(4) 

with uj' = uj + pCl, d = p + q, and 
Cd{a,b,a, 

p 

According to Eq. (|4]) we will observe coincidences be- 
tween Alice and Bob's photons only if the frequency bins 
LUA.B in which they are detected are separated by integer 
multiples of fl: loa + — 2u;q = df2, d € Z. 

The rate at which Alice and Bob will detect photons at 
angular frequencies uja = ujq +i-u' and ljb = — w' -I- dVl 
is proportional to 

Q{u}o+uj\ujo - uj' + dfl\a,b,a,(3) = |cd(a, 6, a, . (6) 

(For a derivation of the proportionality factor, see the 
appendix.) 

Because of the symmetries of Eqs (|1I2I3P the quantity 
Q depends only on the absolute value of the index d (but 
not on the sign of d, nor on luq and w') and on the phase 
difference A = a — /3 (but not on a + /3): 

Q{u}o + u}',uja-uj' + dn\a,b,a,/3) = Q{d\a,b,A) (7) 

= Q{-d\a,b,A). 

The quantities Q{d\a,b, A) obey the normalization con- 
dition 

^ Q{d\a, b,A) = Y, Ma, 6, a, = 1 (8) 

d d 

and satisfy 

0((i = 0|a,a,7r) = 1 , Q(d 7^ 0|a, a, tt) = . (9) 

When the phase modulators are turned off (a = = 0), 
the correlations are trivial, and we have 

Q{d = 0\a = b = 0) = 1 , Q{d^O\a = b = 0) =0. 

(10) 



In the experiment reported below, we chose w' = 0. 
Because the filters Fa,b have a finite bandwidth, the ac- 
tual values of uj' belong to a small interval [— e/2,e/2] 
of width approximately 3 GHz centered on uj' ~ 0. To 
resolve the frequency bins we need both that e < Q and 
that the filter transmission T drops very steeply beyond 
e/2. These conditions are experimentally realized thanks 
to the properties of the filters, see description above and 
Fig. [1] insets (a) and (b). The value of d is chosen by ad- 
justing the reflected frequency of Bob's filter to ujq + dil, 
while the reflected frequency of Alice's fllter is kept flxed 
on ujQ. In Figs [Hand [31 we compare the predictions of Eqs 
(|6I7[) to our experimental results, for d = 0, 1, 2, 3, 4, 5. 

Our experimental estimate, denoted Q, of the quantity 
Q{d\a, b, A) given by Eqs ()6l7p is obtained from the data 
recorded by the time to digital converter (see Fig. [1]) by 
taking the total number of coincidences and subtracting 
the accidental coincidences, and then normalizing by the 
same quantity when d = and a = b = 0: 

Q{d\a,b,A) = 

NMa.b, A)-7V„,(d|a,6,A) 
N^{d = 0\a = b = 0)- Nacid = 0\a^b = 0)'^ ' 

This ensures that Q has the same normalization as Q, see 
Eqs ([HI [9]) , since the coincidence rate is maximum when 
a = b = and d ~ 0. 

In Figs[2]and[3]thc acquisition time per measured point 
was constant, corresponding to a number of coincidences 
approximately equal to lO^P. Experimental values are 
plotted with statistical vertical error bars, which is the 
main source of uncertainty. Note that when d = 5 the 
filter Fb was at the limit of its tuning range, and it may 
have not been perfectly centered on + 5 51, in which 
case there would be a systematic underestimate of Q{d ^ 
5|a,6,A). 

In Fig. [21 the normalized coincidence rate Q{d\a, b, A) 
is plotted as a function of the modulation amplitudes 
a = b (taken to be equal) when A = 0. The number of 
frequency bins that interfere together is approximately 
given by the number of values of d for which Q takes a sig- 
nificant value, and increases when a, b increase. In our ex- 
periment we were able to scan the values a,b G {0, 2.74}. 
When a « 6 « 2.74 there are contributions from 0? = to 
d ~ 5 (and by the symmetry of Eq. ([7|) there should also 
be contributions from c? = — 1 to c? = —5). This shows 
that at least 11 frequency bins are coherently addressed 
by the phase modulators. 

Fig. [3lis obtained by scanning the phase A when a ~ 
b sa 2.74. Note that when A — tt only d — contributes, 
as predicted by Eq. ([9|). The quantity Q{d — 0\a = 
b = 2.74, A) vanishes for specific values of A = A*, 
see the theoretical curve in Fig. [31 This allows us to 
estimate the visibility of interferences through the usual 
formula V = (Qmax - Qmin) / (Qmax + Qmin), where we 
take Qmax = Pid = 0|a = & = 2.74, tt) and Qmin = 
Q{d = 0\a = b = 2.74, A*). From the data reported in 
Fig. [31 we estimate that the visibility is approximately 
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equal to 98%. 



III. NON LOCAL CORRELATIONS 

Now we show how frequency bin entangled photons 
addressed locally by phase modulators can in principle be 
used to realize quantum non locality experiments. That 
is, we want to show that the correlations between Alice 
and Bob's detectors should not be explainable by a Local 
Hidden Variable (LHV) model. 

Our starting point will be the Clauser-Horne [s^l in- 
equality which must be satisfied by any LHV model: 

P(00|AiBi) + P(00|AiB2) + P(00|A2Si) 

~P{m\A2B2) < P{0\Ai) + P(0|Si) , (12) 

where Ai,A2 are two possible settings of Alice's mea- 
surement apparatus and Bi , B2 are two possible settings 
of Bob's measurement apparatus. Each measurement has 
two possible outcomes, and denotes one of the outcomes 
of the measurements. The other outcome does not inter- 
vene explicitly in the inequality. Thus P {00\ AB) denotes 
the probability of both Alice and Bob finding outcomes 
given settings A and B; and P(0|^) = ^2x^(0^1^^) 
denotes the probability of Alice finding outcome 0. Since 
Alice and Bob's measurements are independent, P(0|A) 
does not depend on B. 

In our case the measurement settings will correspond 
to a choice of amplitude and phase applied to the phase 
modulators: Ai = (ai,ai), A2 = (02,02), Bi = (61, 
B2 = {h2,P2)- We will take the outcome in Eq. (fT2|) 
to correspond to the photon being registered in the fre- 
quency bin centered on wq. 

The probabilities P{QQ\AiBi) can be estimated directly 
with our setup as they will be proportional to the num- 
ber of coincidences if both filters Fa.b are centered on 
frequency wq, see Eq. ([TT]) . 

The quantities P(0|^i) and P{Q\Bi) cannot be mea- 
sured directly with our setup, as we only had two detec- 
tors, one on Alice's side and one on Bob's side. However 
we can estimate these quantities by making the follow- 
ing assumption (identical in spirit to the one made by 
Clauser-Horne in [s^ ) : 

P(0|^i) ^ P(0|Bi) = P(00|a = = 0) . (13) 

That is, we assume that the number of photons detected 
by Alice (Bob) in the frequency bin centered on loq when 
Alice (Bob) detector has setting Ai (Bi) is identical to 
the number of coincidences in frequency bins cjq when the 
phase modulators are turned off {a = b = 0). Quantum 
mechanics predicts that this inequality is obeyed, since 
it follows from the symmetries of the correlations Eq. 
(O. the normalization Eqs (j8l9p . and the fact that the 
correlations are trivial when a = b = 0, Eq. (jlO[) . Our 
assumption is that the LHV model also obeys Eq. ([T^ . 



We then insert Eq. (|T3|) into Eq. (|T2)) and divide by 
P(00|a = 6 = 0) to obtain 

P(00|AiBi) P(00\AiB2) . 
P(00|a=6=0) ~^ P{OQ\a=b=a)' 

P(00|A2Bl) P{00\A2B2) ^ r, /-..N 

P(00|a=b=0) P(00|a=6=0) — ' V-^^^ 

If we express that the probabilities should be propor- 
tional to the number of coincidences minus the number 
of accidental coincidences, we obtain the inequality 

S = Q{ujo,uJo\AiBi) + Q{uJo,ujo\AiB2) + 

Q{iOo,uJo\A2Bi) ~ Q{ujn,ujn\A2B2) < 2, (15) 

where Q is given by Eq. (fTTj) . 

On the other hand, a maximally entangled state could 
lead to values as high as 2\/2. In order to investi- 
gate whether Eq. (jl5|) can be violated experimentally, 
we took for simplicity the modulation amplitudes to be 
equal, ai^2 = and numerically optimize the phases 
(^1,2, /3i,2- We considered 4 different values of modulation 
amplitude: 

• ai^2 = bi,2 ~ 0.51 , for which the optimal phases 
are ai = 0, ^2 = 1-42, /3i = 3.85, /J2 = 2.43 . 

• ai,2 = ^'1,2 = 1-01 , for which the optimal phases 
are ai = 0, a2 = 1-02, /3i = 3.65, P2 = 2.63 . 

• ai_2 = ^1.2 = 1.50, for which the optimal phases 
are ai = 0, ^2 = 0.72, /3i = 3.50, P2 2.78 . 

• ai,2 = &i,2 ~ 1.95, for which the optimal phases 
are ai = 0, ^2 = 0.56, /3i = 3.42, ^2 = 2.86 . 

As one can see in Fig. 21 these choices lead to a strong 
violation of the bound 2. 

The observed correlations do not however reach the 
theoretical optima for large values oi a = b. The reason 
why it is difficult to reach the optimal value for large 
values of a = 6 is that the curve P{d ~ 0\a = b. A) is 
more strongly peaked around A = tt for large values of 
a = b; and for the optimal values, A = a — /3 lies on the 
slopes of this peak. 

We have estimated what could be the effect of slight 
errors on a and b (in particular letting b be slightly differ- 
ent of a), and of slight errors in the phases ai,(3i. In the 
experiment the error on /3i was probably larger than the 
error on ai because the phase shifter used to choose 13 
was of lesser quality. Letting a, b vary by a ■ 10~^ around 
the estimated value, a vary by 5 • 10~^ around the ideal 
value, and /? vary by 10 • 10~^ around the ideal value, 
and taking the worst case, we would obtain the curve in- 
dicated in dashed in Fig. |4l These estimates of the errors 
thus provide a possible explanation for the discrepancy 
between the theoretical optima and the observed viola- 
tion of Eq. (fT5|) . 
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IV. CONCLUSION 

In summary we have demonstrated how, using phase 
modulators and narrow band filters, one can accurately 
address in the frequency domain photons belonging to a 
high dimensional entangled state. In view of the proven 
success of sideband encoding for long distance QKD [Tol - 
[22j |. this seems a promising technique for quantum com- 
munication. On the other hand, the class of unitary 
transformations explored in this work is somewhat lim- 
ited, as it depends only on 2 parameters, see Eqs (HIE]). 
However the use of non sinusoidal voltages would give 
rise to much more general families of unitary transfor- 
mations. We hope to explore this in future work. We 
will also focus on studying other Bell inequalities, such 
as the CHSH [H and CGLMP [s^ inequalities. 

We acknowledge support from the European Union un- 
der project QAP (contract 015848), from the Belgian Sci- 
ence Policy under project IAP-P6/10 (Photonics@be), 
from the French Agence Nationale de la Recherche un- 
der project HQNET and from the Conseil Regional dc 
Franche-Comtc. This work also benefits from the Pro- 
gramme International de Cooperation Scientifique PICS- 
3742 of the French Centre National de la Recherche Sci- 
entifique. 



Appendix A: Normalization 

We take the state \uj) to be normalized to 



= (5(0) 



T 
2^ 



We rewrite Eq. (O as 



f + OO 



du)f{uj)\ujQ + a;)A|i.^o - ^) 



(Al) 



(A2) 



where f{uj) takes into account the finite bandwidth of the 
signal and idler photons. We will consider the case where 
/ is a slowly varying function which can be considered 
approximately constant when lj changes by order f2. 
Using Eq. (jAip we have 



p-\-oo 



dx\fi, 



(A3) 



which means that the photon pairs are produced at the 
rate 



1 /■+°° 



(A4) 



Consider now the operator that projects onto a fre- 
quency bin of width e: 



n„ 



uj — e /2 



duj'luj') (w'l . 



(A5) 



The rate of coincidences in two frequency bins of width 
e symmetrically spaced on either side of ujq is 



T 

= 2^ 



J+e/2 



-e/2 



dx\f{x)\' 



(A6) 



The quantum state after it passes through the phase 
modulators is (see Eq. ^) 



J doj \ujo + oj + pn)A\i^o — uj + gO) B 

Up{a,a)Ug{b,(3)f{w) (A7) 
/ du>' ^ \luq + u>') a\(^o — + dft) B 

d 

Up{a, a)Ud-p{b, P).f{Lo' - pn)^ (A8) 
/ dw7(w') Iwo + ^^OaI^o - + dVl)B 

d 

(^Upia,a)Ud-pib,P)^ (A9) 
/ da;'/(a;') ^ \uJo + uj')a\^o - uj' + d^l)B 



Cd(a, 6,q;,/3) , 



(AlO) 



where in obtaining line (|A9[) we have used the fact that, 
for fixed a, &, UpUd-p decreases rapidly with p and that / 
varies slowly so that /(oj' — pVl) ~ /(i^') for the relevant 
p. The coefficients arc given by Eq. ([S]). 

The rate of coincidences in two frequency bins dis- 
placed one with respect to the other by dfi is 



mK>+^'K.-^'+dnW) 

2 T r^+(-/2 ,1 fr\\2 



\c,[aA^)?^!:::'/,Mn^)\ 

= \cMb,^)?Re{o^')T 



(All) 



Since R^{uj') varies slowly with w', the quantity 
Q{d\a,b, A) ^ | c^; (a, &, A) p can be estimated as described 
in Eq. pT|) . 

Note that the normalization condition Eq. (jS]) 
Y^^\cd{a,b, A)\^ = 1 expresses the fact that the trans- 
formation l^*') is unitary and that no photons are 
lost in the process. 
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FIG. 1: Experimental setup. The quasi-monochromatic pump 
laser (PUMP) creates photon pairs in the Periodically Poled 
Lithium Niobate waveguide (PPLN) and is then removed by 
a drop Filter (F). Alice and Bob's photons are selected pas- 
sively by a 3dB splitter. The photons then pass through 
electro-optic Phase Modulators (PMa.b)- The phase modu- 
lators are driven by a 12.5 GHz Radio Frequency Generator 
(RF GEN) whose output is controlled by variable attenuators 
(Va,b) and phase shifters {a, j3). RF isolators (not shown) 
guarantee the independence of Alice and Bob's settings (75 dB 
isolation). Individual frequency bins are selected by narrow 
band filters which consist of a Fiber Bragg Grating (FBGa.b) 
preceded by a Circulator (Ca,b). The photons are finally de- 
tected by Avalanche Photo-Diodes (APDa.b) and the elec- 
tronic signals sent to a Time to Digital Converter (TDC) to 
perform a coincidence measurement. A typical figure recorded 
by the TDC, consisting of a background due to accidental 
coincidences and a narrow peak when the photons arrive in 
coincidence, is shown. Typical values for the Signal to Noise 
Ratio (SNR) (number of non accidental coincidences divided 
by number of accidental coincidences) are SNR ~ 100. Inset 
(a) shows a typical transmission spectrum of the narrow band 
filters. Inset (b) illustrates how the phase modulators gener- 
ate new frequencies: if a broadband source passes through 
the narrow band filter (as in (a)), and then through a phase 
modulator, one obtains the spectrum of (b). The height of 
the peaks is given by the norm square of the coefficients in 
Eq. ([T]). When a ~ 2.74 (corresponding to the maximum RF 
power produced by the source) one can see eleven frequency 
bins corresponding to the order of the peak p, see Eq. ((TJ, 
ranging from p = —5 to p = -1-5. In insets (a) and (b), A is the 
wavelength and T the transmission. Note that the horizontal 
scale is different in insets (a) and (b). 
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FIG. 2: Theoretical predictions (curves) and experimental 
measurements (with error bars) of the normalized coincidence 
rate Q{d\a, b, A) = Q{d\a, a, 0) when A « and the amplitude 
a is scanned, for d = 0, 1, 2, 3, 4, 5. The experimental measure- 
ments are plotted entirely in terms of measured quantities, 
and do not depend on any adjustable parameters. Values of 
a are deduced from measures of the RF power. Horizontal 
error bars are due to the limited resolution of the power me- 
ter used (we assumed a relative uncertainty on a of 10"^). 
Vertical error bars are statistical. 
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FIG. 3: Theoretical predictions (curves) and experimental 
measurements (with error bars) of the normalized coincidence 
rate Q{d\a, b, A) = Q(d|2.74, 2.74, A) when a b ^ 2.74 and 
the phase A is scanned, for d = 0, 1, 2, 3, 4, 5. To plot the ex- 
perimental measurements we used the value indicated by the 
mechanical phase shifter - call it x -, and converted it to a 
phase value using the relation A = nx + u. Parameters /i and 
1/ were adjusted to get a good fit with the theoretical predic- 
tions. Horizontal error bars are due to the limited resolution 
of the phase shifter used (we assumed an absolute uncertainty 
on A of 5 ■ 10"^-ad). Vertical error bars are statistical. 
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FIG. 4: Violation of the inequality Eq. (|15|) as a function 
of the modulation amplitudes ai,2 = foi,2 = a. The experi- 
mentally determined values of S are given by the left hand 
side of Eq. (|15p . They include statistical error bars (vertical 
axis) and RF amplitude error bars (horizontal axis) . The top 
curve shows the theoretical evolution of the maximal value of 
S when ai,2 and /3i,2 are numerically optimized. The dot- 
ted curve takes into account possible errors on ai,bi, ai, Pi: 
it shows theoretical predictions when phases and amplitudes 
are allowed to fluctuate around optimal values (see text for a 
detailed description). 



